function align_axis
fold=fopen('grid.bdf','r');
fnew=fopen('grid+.bdf','w');
while feof(fold)==0
    entry=fgetl(fold);
    if length(entry)<2
        continue
    end
    x=sscanf(entry(25:32),'%f');
    y=sscanf(entry(33:40),'%f');
    z=sscanf(entry(41:length(entry)),'%f');
    len_x=length(x);
    len_y=length(y);
    len_z=length(z);
    if len_x==2
        entry(25:32)='    0.0 ';
    end
    if len_y==2
        entry(33:40)='    0.0 ';
    elseif y<0.01
        entry(33:40)='    0.0 ';
    end
    if len_z==2
        entry(41:48)='    0.0 ';
    end
    fprintf(fnew,'%s\n',entry);
end
fclose('all');
!del grid.bdf
!rename grid+.bdf grid.bdf
end